Weak disorder expansion for localization lengths of quasi-lD systems 
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A perturbative formula for the lowest Lyapunov exponent of an Anderson model on a strip is 
presented. It is expressed in terms of an energy dependent doubly stochastic matrix, the size of 
which is proportional to the strip width. This matrix and the resulting perturbative expression for 
the Lyapunov exponent are evaluated numerically. Dependence on energy, strip width and disorder 
strength are thoroughly compared with the results obtained by the standard transfer matrix method. 
Good agreement is found for all energies in the band of the free operator and this even for quite 
large values of the disorder strength. 
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The Anderson model describes the generic behavior of 
the motion of an electron in a disordered solid. In the 
one-dimensional (ID) situation, rigorous proofs of strong 
localization have been given [ij. Moreover, the local- 
ization length has been calculated 0, Q for weak disor- 
der and its inverse, the Lyapunov exponent, is given by 
7 « wV(96sin2 k) = w"^ / - 2A E"^) , where E ^2cosk 
is an energy in the Bloch band of the unperturbed op- 
erator (away from band center and edges) and w is the 
disorder strength in the normalization discussed below. 
For the 2D case, scaling theory |J| predicts also strong 
localization with an disorder dependence of the localiza- 
tion length of the form exp(l/w^). This was confirmed 
by high-precision numerical studies based on the transfer- 
matrix method (TMM) 0,0]. In 3D, a transition to the 
so-called weak localization regime with diffusive motion 
is expected for low disorder strength and energies in the 
band of the free operator. A rigorous proof of strong 
localization in 2D and 3D exits only for the band edges 
and at high disorder 0, |1| and, in particular, not for en- 
ergies in the band and small disorder in the 2D situation. 
No rigorous breakthrough results are known for the weak 
localization regime. 

In order to approach the higher dimensional cases, a 
detailed understanding of the quasi-lD situation, i.e. an 
infinite wire with many channels, is of crucial importance. 
The TMM d, Q is here a very reliable tool for a numer- 
ical study of the inverse localization length, namely the 
smallest Lyapunov exponent. Thouless argued that it is 
proportional to uP' jL where L is the number of channels 

. The Dorokhov-Mello-Pereyra-Kumar theory gives a 
microscopic derivation of this behavior (see ^3 for a re- 
view), but is not based on a calculation starting directly 
from an original model and hence does not allow to study 
finer properties such as the energy dependence of the lo- 
calization length. A perturbative analytical calculation 
of the smallest Lyapunov exponent of an Anderson model 
on a strip was recently given by one of the authors . 
The techniques of that work also allow to deal with other 



models having symplectic transfer matrices. 

In this letter, we first recall the perturbative formula 
from [Tl| and then compare it numerically with the 
TMM. Our main result is that the perturbative formula 
works remarkably well for all but a discrete set of energies 
and, quite surprisingly, relatively large values of the dis- 
order strength. This allows to understand the rich struc- 
ture of the energy dependence of the smallest Lyapunov 
exponent. We also study various quantities associated to 
the random dynamical system underlying the TMM. 

Let us begin by recalling that the Anderson Hamil- 
tonian on an infinite strip of finite width L is given by 



H 



(1) 



{x,y> 



with tight-binding states |x) at x = (n, to) where ti G Z 
and TO = 1, . . . , L. The hopping is between nearest neigh- 
bors (x, y) only and with periodic boundary conditions 
in the m-direction. The w(x) G ^] ^"^^ independent 
and identically distributed random variables with (for 
sake of concreteness 0|) uniform distribution so that 
E [^(x)^] = For simplicity, let us also choose L even. 
The Schrodinger equation HvP = E'^ is rewritten in a 
recursive form using the 1L x 2L transfer matrices 



T{n) = 



w'V{n) 
1 
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(2) 



where A/, is the discrete Laplacian in the transverse di- 
rection with periodic boundary conditions and V(n) — 
diag(w(n, 1), . . . , t)(n, i)). The transfer matrix is sym- 
plectic, namely, T(n)*JT(n) = = ^ ^ Asso- 
ciated with this family of random matrices are the Lya- 
punov exponents 71 > . . . > 7l > defined via 
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for p = 1, . . . ,L, where denotes the p-fold exterior 
product. They are self- averaging quantities so that an 
average over the disorder configurations may be taken 
before the large N limit [l^ . 

The first aim is to bring the transfer matrix at w = to 
its symplectic normal form. Let /i; = —2 cos {27tI / L) — E , 
I = 0, . . . , L — 1, denote the eigenvalues of — E. 
Clearly, /u/ — ^l-i so that all eigenvalues /z/ except for 
I — 0, L/2 are doubly degenerate. If |^/| < 2, we call 
it an elliptic eigenvalue and define its rotation phase rji 



by m = 



If on the other hand, 



> 2, 



we call it hyperbolic and define its dilation exponent r/i 
by = e 



e '" . There are — 1 energies E within 
the band [—4, 4] of the free strip Laplacian for which the 
parabolic case — 2 occurs. These energies corre- 
sponding to interior band edges are for now excluded, 
but will be further discussed below. Using the corre- 
sponding eigenvectors of Al it is possible to construct a 
symplectic matrix M such that 



M-iT(n)M = R[l + wP(n) 



(4) 



The symplectic matrix R is built from elliptic and hyper- 
bolic rotation matrices Re(?7), Hhirj), respectively given 
by 

/ COST] — sin?7\ / cosh 77 sinh?] \ 
I sin 77 cos?7 I ' I sinhjy cosh 77 J 



(5) 



More precisely, the 4 entries of R at {1,1), {l,L + I), 
{I + L,l) and {I + L,l + L) form the matrices Re(?7;), 
R/i(77i) depending on whether fii is elliptic or hyper- 
bolic. All other entries of R vanish. The matrix P(n) = 

M^^ ^ V(n) ) nilpotent and in the Lie algebra 

of the symplectic group. 

These free modes naturally group themselves into 
channels which are the even-dimensional subspaces of 
M."^^ rotating under R with the same frequency. In our 
situation, there are i/2 -I- 1 such channels indexed by 
I — 0, . . . , L/2. The Zth channel is given by the compo- 
nents I, L — I, L + I, 2L — I oi R^^ and we denote the 
corresponding projection by tt;. The channels I = 0, L/2 
are simple, while all others are doubly degenerate. 

Let us introduce a symplectic frame u = (ui, . . . 
to be a set of orthonormal vectors in R^^ satisfying skew- 
orthogonality {up\J\uq) = for all p,q — I, . . . , L. Given 
an initial symplectic frame u(0), new frames u{n) are 
constructed iteratively as follows: apply R [1 -I- wP{n)] 
to each vector of u{n — l) and then use Gram-Schmidt or- 
thonormalization procedure in order to obtain u{n). This 
gives a random dynamical system on the space of sym- 
plectic frames (which is isomorphic to the L-dimensional 
unitary group). The advantage of the basis change Q 
is that the discrete-time dynamics of frames at w = is 
simply given by rotations. A weak random potential per- 
turbs this simple dynamics in an analytically controllable 
way. 

Important in our perturbative formula will be the 
weight Pp^k{n) = {up{n)\-Kk\up{n)) of the pth frame vector 
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FIG. 1: Lyapunov exponent 7l as function of disorder 
strength for energy E — 1.1 and strip widths L — 20, 50, 100. 
Lines have been computed from ||7J with 0.05% error, symbols 
denote results of TMM calculations (with 0.5% error if there 
is no error bar). Error bars are drawn on one side only. 



in the kth channel at iteration n, as well as its Birkhoff 
mean 



lim 



^7 H Pp,ki^ 



(6) 



n=0 



Each matrix Pp^k{n) is doubly stochastic, namely, the 
sum over the channel index fc = 0, . . . , L/2 equals 1 while 
the sum over the frame vector index p — 1 , . . . , L is equal 
to the degeneracy of the fcth channel. The latter fact 
is related to the symplectic structure of the frame and 
[7ri,J] =0. In a similar fashion, one can define other 
Birkhoff averages such as {pL,iPL,k)- 

The rigorous perturbative formula for the smallest 
Lyapunov exponent within the band is then, under a 
hypothesis on the incommensurability of the rotation 
phases which excludes energies with Kappus-Wegner- 
type anomalies jl^ and interior band edges, given by 



96 L ^ sin ?/; sin rjk 



{PL,1 PL,k) + Oiw^) , (7) 



where the sum runs over elliptic channels only (actually, 
for hyperbolic channels k, pL,k{n) almost vanishes for n 
large enough as will be discussed below). For a single 
channel, this expression reduces to the perturbative ID 
result given above if one sets Ai = in (O and then 
Po = E so that r]o = k. We emphasize that only the av- 
eraged channel weights of the last frame vector ul enter 
expression Q. We also remark that the dependence of 
the error term on L and E remains unspecified. 

For a numerical study of Q), one first evaluates the 
Birkhoff averages {pL,iPL,k) as in © by generating ran- 
dom transfer matrices just as in the TMM 15] . Report- 
ing this into (0 and neglecting the 0{w^) term gives the 
perturbative values plotted in Fig. ^ (as well as Fig. [21 
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FIG. 2: Lyapunov exponent jl as function of energy for dis- 
order w — 0.2, 0.5, 1.0 and strip width L — 20. Solid lines have 
been computed from ||7|l with 0.05% error, symbols denote re- 
sults of TMM calculations with 0.5% error. For w = 0.5 the 
dashed solid line was calculated from Eq. ||HJ| and for w — 0.2 
from Eq. ©. 



FIG. 3: Log-log plot of the reduced localization length 1/I/7l 
as function of inverse strip width 1/L for energy _B = 1.1 and 
disorders w — 0.2 (o), w = 0.5 (□), and w = 1.0 (o). Full 
symbols denote results of TMM calculations with 0.5% error, 
open symbols have been computed using l|7| with 0.05% error. 
The lines serve as guide to the eye. 



andOlbelow). This was done for a typical energy i? = 1.1 
away from internal band edges such that ^ 2 for all 
I. For comparison, we also plot the value of 7^ as evalu- 
ated by the standard TMM. Very good agreement is ob- 
tained even for rather large values of w. For very small 
w, the numerical convergence of the TMM estimates for 
7l is computationally intensive whereas the convergence 
of the average channel weights needs about a factor 10'^ 
less iterations. More stable results can hence be obtained 
at a fraction of the computational cost. For L = 20 and 
E — 1.1, the validity of the perturbative formula breaks 
down at about w « 3.0, for large L a bit earlier. The 
breakdown happens in the region of crossover from quasi- 
ID to 2D behavior, i.e. Lji^ « 1. 

Fig- HI shows the energy dependence of 7^ and its per- 
turbative approximation Q . The oscillatory behavior is 
very well reproduced by the internal band edges corre- 
sponding to some eigenmode of A^,, that is the peaks 
lie precisely on energies E where — 2 for some I so 
that sin 77; — 0. As the validity of ^ in their vicinity is 
restricted, the singularities of the perturbative approxi- 
mation are artificial. In fact, one has to adapt the anal- 
ysis of ^5 using the normal form of |T3| to this higher 
dimensional case. This will be done elsewhere. 

Fig.^shows the strip width dependence of the reduced 
localization length in a plot usually used to infer the 
universal 2D scaling function. The scaling hypothesis 
^ J8J states that there is a unique scaling function F 
such that one can find a function ^{w), also called the 
2D localization length, so that all the data 7l(u') satis- 



fies log (^) = F 



Hence Fig. O allows 



log(^ 

to read off a part of F which is usually difficult to de- 
termine numerically because the TMM converges badly 
for very small w. Note that all the data of Fig.|3|lie well 
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FIG. 4: A sample of numerical results for the average channel 
weights {pp,k) as function of the channel number k for various 
frame vectors (index p). The strip width is L = 50, energy 
E = 1.1 and disorder w = 0.5. HyperboHc frame vectors with 
p < 17 have a large {pp,k) ~ 1 for a specific k, whereas elliptic 
frame vectors p > 18 all have nearly a uniform distribution on 
elliptic channel indices 9 < A: < 25. The lines serve as guide 
to the eye. 



inside the quasi-lD regime where Ljl 1. 

We now analyze further quantities of the TMM random 
dynamical system. Fig. 01 shows typical numerical results 
on the average channel weights. The first frame vectors 
align with the hyperbolic channels. More precisely, they 
fill them one after another in order of decreasing dilation 
exponents. Thus ui completely aligns with the expand- 
ing direction of the most hyperbolic channel correspond- 
ing to the fundamental of A^. Vectors p = 10, 11 fill the 
doubly degenerate sixth channel. Due to symplectic and 
orthogonal blocking, the remaining frame vectors have 
to be in the elliptic channels. In fact, they have a more 
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FIG. 5: Contour plot of the numerically evaluated correla- 
tion function {pL,ipL,k) — {pL,i}{pL,k) of the Lth frame vector 
plotted as function of channel indices I, k for energy E — 1.1, 
disorder w — 0.25, and strip width L = 50. The correlation 
for hyperbolic channel numbers Z, fc < 8 is zero. 

or less uniform distribution over the elliptic channels (in 
particular, also the last one p = 50), only the first elliptic 
frame vectors p = 18, 19 give more weight to channels 
close to the internal band edges. Moreover, further nu- 
merics shovir that all these distributions are nearly inde- 
pendent of the disorder strength for w sufficiently small 
enough. Indeed, one can set up a BGGKY-type hierarchy 
for these distributions • The fact that the distribution 
over the elliptic channels is nearly uniform is a result of 
the mixing properties of the random potential, because 
the free dynamics R gives no preference to any of the el- 



liptic channels. Fig. [S] shows that the correlations in the 
calculation of Birkhoff averaged channel weights are very 
small. Hence, as further approximation, one may factor 
the stochastic matrix in Q. 

Resuming, the numerical results imply that the dis- 
tribution of the last frame vector is uniform on the el- 
liptic channels, the correlations are weak and all this is 
uniformly for small w. Hence a good approximation for 
the coefficient in {TJ should be obtained upon replacing 
{pL,i PL.k) by (Le)~^ if Le is the number of elliptic chan- 
nels. This gives 

w2 1 ^ 2 - 5i^k 
95 L L-^ sm rji sm rjk 

where the sum still runs over elliptic channels only. This 
fits well with the results of Fig. ^ and gives a good ap- 
proximation as shown in Fig.|21for w = 0.5, albeit not as 
good as Eq. {Tj). For large i, one may furthermore ne- 
glect the 5i^k in ® approximate the discrete sum by 
a Riemann integral. Hence, we infer for E > Q and large 
L (but not too large so that E stays away from internal 
band edges) 

12-^ V"'''^^''^^ ~ V4- (2cos77 + S)2 j 

(9) 

where we set r\E = arccos(l — E/2). Indeed this elliptic 
integral (of first kind) can be evaluated numerically. Of 
course, this does not give the rich oscillatory structure 
for finite L anymore as demonstrated in Fig. [3 for the 
case w — 0.2. 
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